Van-Hove singularities effects on the upper critical field 
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We present a study of the superconducting pairing susceptibility Kr(r) for a two-dimensional 
isotropic system with a strong power-law divergence in the density of states N(e) ~ e _1+ ^, b > 1. 
We show that the pair propagator has the scaling form, Kt(t) = r b ~ 3 f 1 (T 1/ ' ! Y). An anomalous 
short-range behavior is found, leading straightforwardly to positive curvature in the upper critical 
field, for 6 < 2 and to a zero temperature divergence, H C 2 ~ T~ 2+ b, for b > 2. 
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PACS numbers: 74.60. Ec 

Photoemission experiments in Copper-Oxide super- 
conductors [Q have provided direct evidence for the ex- 
istence of an extended saddle point in the Cu0 2 plane 
bands and consequently, a strong divergence in the den- 
sity of states, N(e) ~ (e — e v h)~ a , for energies close to 
e v h- Some authors have claimed that the high supercon- 
ducting critical temperatures of the cuprates could be 
explained taking into account this divergence 

One of the most surprising properties of these mate- 
rials is the upper critical field, which has been obtained 
in magnetoresistance experiments down to very low tem- 
peratures in the case of overdoped Tl 2 Ba 2 Cu06+5 
Bi 2 Sr 2 CuO y H and underdoped YBa 2 Cu307_i 
very unusual if c2 (T) curve is observed, with very strong 
positive curvature and no evidence of saturation at low 
temperatures. This behaviour contrasts strongly with 
the weak coupling BCS result @ which predicts an ap- 
proximately parabolic shape for the 7? c2 curve. 

Recently, Abrikosov has proposed (|] that these 
anomalous _ff c2 curves reflect a dimensional crossover to 
quasi-one dimensional superconductivity due to the pres- 
ence of flat regions in the energy dispersion, that is, ex- 
tended saddle points. In Abrikosov's approach, the two 
extended saddle points in the energy dispersion are re- 
placed by two one-dimensional linear energy dispersions 
eife) = vi-Qx and e^fay) = u 2 -Qy This model is equiva- 
lent to a system of two transverse chains and in this case, 
the density of states loses completely its strong energy de- 
pendence. Furthermore, it is not surprising that he finds 
a dimensional crossover in i? c2 . In this paper, we ar- 
gue that these curves reflect not a dimensional crossover, 
but the strong energy dependence of the density of states 
which results from the presence of these extended saddle 
points. In the following, we present a study of the su- 
perconducting pairing susceptibility for a isotropic two- 
dimensional system with a strong power-law divergence 
in the density of states. As shown by Gorkov JjJ, this 
two-particle correlation function determines the shape of 
the superconducting transition i/ c2 (T) of a type- II su- 
perconductor. 

The superconducting transition is characterized by 
the vanishing of the gap function A(r,r'), defined as 



A(r,r') = V(r — r')(ipi(r)ipi(r')). In the particular case 
of a local pairing interaction, V(r — r') — g S(r — r'), 
we obtain as usual the s-wave gap function, A(r,r') = 
A(r)6(r — r'). In the following, /j, = c= e = fcs = l- 
In the vicinity of the superconducting transition curve, 
the gap parameter is small and a perturbation expansion 
in powers of A leads to the semi-classical linearized gap 
equation RF 



A(r) 
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dr'K (r' - r)e 42A(rHr '~ r) A(r'). (1) 



where Kp{r) is the fermion pair propagator in real space 
for a given temperature T = 1/(3, in the absence of the 
external field and the pairing interaction g and is defined 

as 



Kp(r',r) 



G- U {r' ,r)Q u (r' ,r) 



(2) 



where the Matsubara Green's function describes the 
normal state in the absence of magnetic field. Using 
Kramers-Kronig relations, Kp{r) can be rewritten as 



Kp{r) = - J duitaah(/3u)/2)A(r,u)B(r,-uj), 



(3) 



with A(k,w) = lmG R (k,uj) and B(k,u) = ReG R (k,uj) 
where G R (k,ui) is the retarded Green's function in the 
absence of magnetic field and pairing potential. A(r,oj) 
and B(r,oj) are the respective Fourier transforms. A 
non-local V(i — r') may lead to a d-wave gap solution and 
a slightly modified gap equation. One can show that the 
upper critical field probes the behaviour of the Cooper 
pair center of mass and the internal symmetry of the gap 
function is irrelevant as long as the thermal and magnetic 
lengths are much larger than the interaction range. 

In a bidimensional system, a van-Hove singularity 
(VHS) in the density of states results usually from the 
presence of a saddle point in the energy dispersion e(k). 
In the case of an extended saddle point, e(q) ~ q" — q™, 
where q = k — k v h , this leads to a power-law divergence 
in the density of states N(e) ~ e~ 1+ ~ + ~ . Such form for 
the extended saddle point is not only indicated by the 



1 



direct probing of the energy dispersion using the angle 
resolved photoemission technique, jjj but also by numer- 
ical work on the Hubbard model. For instance, Quantum 
Monte Carlo work by Imada and collaborators [jl0| in the 
Hubbard model has found an extended saddle point with 
a quartic q y dependence at (0, 7r). The effect of such a 
divergence on the pairing susceptibility and the super- 
conducting phase diagram is the subject of this paper. 
Clearly, a system with a saddle point is not isotropic. 
However, in order to simplify the problem, we adopt the 
isotropic dispersion relation: 



e(fe) - e vh = a ■ sign(g)|q| t 



(4) 



where q = k — k v h- The influence of anisotropy in the 
semi-classical upper critical field is well studied. [ pd|Jl2| 
Anisotropic two-dimensional systems have tipically open 
warped Fermi surfaces or elliptical closed Fermi surfaces 
in the case of small particle number. For a system with 
an elliptical Fermi surface, in the case of a transverse 
magnetic field, it is simple to show that the normalized 
upper critical field follows the parabolic-like BCS curve. 
|f3| For an open warped Fermi surface, the behavior of 
the H C 2 curve is determined by the relation between T co 
and the small t y modulation of the Fermi surface. [ ^3|Jl2| 
If t y <C T co , H C 2 will diverge at a finite temperature, 
reflecting a reduction of the effective dimension of our 
system induced by the magnetic field. However, 
if t y 3> T co , a BCS-like parabolic H c2 curve is obtained. 
A reduction of t y enhances the zero temperature critical 
field, H co relatively to the zero field critical temperature, 
T co , but, as long as the relation is valid, the reduced up- 
per critical field (H C 2/H co as a function of T/T co ) remains 
unchanged. One may therefore conclude that, unless a di- 
mensional crossover is present, the reduced upper critical 
field shows very little sensitivity to anisotropy. Another 
important point is that the contribution to superconduc- 
tivity of the extended saddle-point region is much larger 
than the contribution of the other regions of the Fermi 
surface, which can therefore be neglected. These facts 
motivate the choice of an isotropic model for our study. 

Note that, for quasi-2D systems and magnetic field ap- 
plied along the planes, the dimensional crossover is from 
3D to 2D. In the case of transverse fields, the crossover 
is from quasi-2D to quasi-lD superconductivity. One 
should note, however, the H C 2 divergence results from a 
mean-field analysis and fluctuations modify this behav- 
ior greatly in the last case. In fact, saturation should 
arise at low temperatures due to fluctuations, reflecting 
the well known impossibility of a superconducting state 
in one dimension. Furthermore, as recently shown by 
Lebed et al [|l2|, saturation should also be observed due 
to Pauli pair breaking. Therefore, it seems unlikely that 
a dimensional crossover could explain results obtained by 
Mackenzie and others 0g. 

We assume that the VHS is pinned at the Fermi level, 
that is, k F = k v h- We will comment on the pinning 



assumption at the end of the paper. The density of states 
for the above model is N(e) ~ d~ x ' b b~^{e — e v h) Let 
us assume for now that b is an odd integer. 

For this simple model, we can compute the spectral 
function 




2k F 
nr 



sign(w) +kf 



(5) 



and the retarded Greens function G R (r,ui), since 
G R (q,u>) is a meromorphic function in the complex q- 
plane. Note that A(r,u>) = ImG R (r,ui) and B(r,u) — 
Re G R (r,uj). After some lengthy but straightforward al- 
gebra, one obtains the following expression for the pair 
propagator, 



K (r) 



r b-3p 
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with 



F[X] 



2k F 1 f°° tanh[(o,A) b ] f 1 . 
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When X > 1, F[X] - X b - 2 and for X < 1, the func- 
tion is exponentially small. Note that the thermal length 
is given by £t ~ (a/T) 1 / 6 . The pair propagator for dis- 
tances smaller that the thermal length is approximately 
given by Kx{r) ~ r~ 1 T _1+ ^ and therefore, it diverges 
as the temperature goes to zero. We will show that this 
will lead to a zero temperature divergence in the upper 
critical field. Note that no Debye-like frequency cutoff 
was introduced in the previous integral. This procedure 
is valid as long as the temperature provides a smaller cut- 
off in the integrand, that is, T 1 / 6 <C w c . This reflects the 
well known reduction of the isotope effect in the van-Hove 
scenario. jl5| 

The zero field critical temperature is obtained from the 
equation 



1/5 



k F 1 



7T aM h b J 



du) tanh(/3w/2)w b 



i-2 



which leads to 



k F a b 
~b- 1 



i i &/(*>-!) 



(8) 



(9) 



This result and the role of the frequency cutoff can 
be qualitatively understood using the usual BCS rela- 
tion for the critical temperature T co ~ £j c e -1 ' ( iV ( e )) T c<.fl ) 
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FIG. 1. H C 2 curves for 6=3, 5 and 7, obtained numerically 
from Eq. [h]. Inset: The same curves in a logarithmic scale, 
showing clearly the low temperature power-law behavior. 



where (N(e))T ao is the thermal averaged density of states 
(N(e)) T ~ J" de(df/de)N(e) and / is the Fermi distribu- 
tion function. In this case, (N(e)) Tco ~ a~ 1 / b b~ 1 T^ a ~ 1 

and therefore ln.(oj c /T co )T c b ~ l/g. In the weak cou- 
pling limit, one can neglect the logarithmic correction 
and thus, the above dependence for the critical temper- 
ature is reproduced. In the usual case of an extended 
saddle point, this broadening argument leads correctly 
to a transition temperature T co oc g 2 . The en- 

hancement of the critical temperature is clearly bounded 
by the equivalent of the Debye Temperature in this prob- 
lem, that is, Tu m ~ uj£), where lod is our cutoff in fre- 
quency. The energy dispersion as given by Eq. ^ may also 
be limited to an energy range lo c < ujd and in that case 
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This dependence on the extent 



of the anomalous energy dispersion could offer an expla- 
nation for the low critical temperatures of, for example, 
Bi2201 1 1 S | and Sr2Ru04_| 19|_. Photoemission experi- 



ments in these materials |1S ,|19[ have found a VHS near 
the Fermi level, but also a smaller extent of the flat re- 
gions in the energy dispersion. A similar thermal broad- 
ening argument can be applied to the zero temperature 
slope of the H c2 curve. 

The analytical determination of the upper critical 
curve for the complete temperature range is a difficult 
task. So, we obtain the H C 2 curves by numerical solution 
of the gap equation and study its behavior analytically 
only at low temperatures or close to T c . For numeri- 
cal purposes, it is more convenient to work with the gap 
function in a mixed representation. If one chooses the 
Landau gauge A = (0, Hx, 0) and makes use of the de- 
generacy of the gap function, one can rewrite Eq. [I] as 

A(x) = g J dx'K p [x' - x, -H(x + x')]A(x'), (10) 
where A(x) is the y integrated gap function and 



Kp(x,k v ) is the Fourier transform of Kp(x,y). At zero 
temperature, the gap equation simplifies to 



A(r) ~ g dr' 



«<*(-?=, -£=) 



A(r') 



(11) 



where <f> is the magnetic phase acquired by the Cooper 
pair, which is independent of the magnetic field if r 
is written in magnetic length units. Note that this 
form for the gap equation is independent of our gauge 
choice. The numerical gap solutions show a perfect scal- 
ing A(x) — F(x/VJl), that is, all gap solutions fall into 
an universal Gaussian curve (see Fig. ^|), if the x-axis 
unit is the magnetic length and therefore, with the vari- 
able change x — x/\TH in the previous equation, the 
gap function becomes independent of H and we obtain 
the low temperature scaling of the upper critical field, 
H c2 (T)~T-*+i. 

In Fig. 0, H C 2 curves for several values of b, obtained 
numerically from Eq. ^ are displayed. These curves are 
characterized by a strong divergence of the upper critical 
field as T — > and linear behavior close to T co . In the 
inset, the low temperature scaling is clearly observed in 
a log-log scale. This behavior is clearly distinct from 
a dimensional crossover in H C 2 which would lead to a 
divergence even in a log-log scale. 

While Eq. |(j for the pair propagator was derived for 
odd integer b, we believe this equation is qualitatively 
correct for any value of b > 1. Clearly, Eq. || for the 
spectral function is valid for any b and one can show 
that the pair propagator for b > 1 will have the same 
qualitative short and long range behavior as that given 
by Eq. [| For 6 = 1, with the introduction of a cut- 
off, we recover the usual BCS results and, in particular, 
H co ~ Tl Q . For 1 < b < 2, F[X] ~ const, if X < 1 
and therefore, the pair propagator shows a different short 
range dependence, Kt{t) ~ r h ~ 3 . The pair propagator 
does not diverge as we decrease the temperature, and 
with a scaling argument, we can show that now the zero 
temperature critical field is finite, l/g ~ Hil b ^ 2 and 
thus, H co ~ Tco b '. 

The low temperature dependence of H C 2 can be ob- 
tained expanding the pair propagator in powers of T, 



[K T (r) - K (r)]/r 



6-3 



_( rT i/!y 



(121 



and following Gorkov M , one obtains H C 2 (T) — H C 2 (0) ~ 
— r 2c / h . Curiously, a power-law low temperature depen- 
dence of H C 2 has also been suggested by Kotliar and 
Varma [pfjf as a consequence of a zero temperature criti- 
cal point. This dependence, in our picture, results from 
the scaling form of the pair propagator as given by Eq. |(], 
but the value of c depends on the specific form of the in- 
tegrand of Eq. One knows that when b — > 1, the usual 
expression for the pair propagator should be recovered, 
which is the one given by Eqs. [3] and [f] only with the 
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FIG. 2. The numerically obtained gap solution at different 
points of the transition curve shown in the inset (6 = 5). All 
solutions are Gaussian. 



sine function in the integrand |9|,|2l|. For b > 2, the 
exponential term in Eq. [7| dominates and the sine contri- 
bution becomes irrelevant. Therefore, when 6 — > 1, the 
low temperature behavior should be determined by the 
sine term and as b goes away from 1, the exponential 
term should take over. With this assumption, c can be 
determined and the result is c = (2 — b)/2, when b ^ 2 
and c = (3 — b) /2, when b — ► 1. 

The results obtained up to now can be collected 
into a equation similar to the usual one, plf 1/g — 
J drKp{r)e~ r H with a qualitative pair susceptibility 
given by 



K T {r) 



1 



(rTV6) 



r3 b sinh[(rTi/ft) d 



(13) 



with c = and d = b - 2, if b > 2. If 1 < b < 2, c = d 
with c having the behavior described above in order to 
reproduce the low temperature dependence of the upper 
critical field. In particular, c = 1, if b = 1 and the usual 
BCS equation is recovered If b — > 2, c — > 0. In 

Fig. ||, i? C 2 curves obtained with this qualitative kernel 
are displayed. A drastic transformation from conven- 
tional parabolic- like curves (obtained with c = (3 — 6) /2) 
to curves with strong positive curvature (obtained with 
c = (2 — b)/2) is observed as the low temperature expo- 
nent 2c jb goes from 2 to 0. 

In Fig. |3J, the experimental H C 2 points obtained by 
Mackenzie et al for Tl 2 Ba 2 Cu06+,5 @] are also displayed 
and fitted with our qualitative i? c2 curves. Note that this 
is a one-parameter fit (c — (2 — b)/2), since the normal- 
ized curves do not depend on the coupling constant g. An 
impressive agreement is observed for 2c/6 = .45, which 
according to the picture presented in this paper, implies 
that the density of states diverges as N(e) ~ e~- 28 . Most 
photoemission experiments have found saddle points with 
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FIG. 3. The qualitative normalized H C 2 curves for 
1 < b < 2. As the dispersion relation changes from linear 
to quadratic, the upper critical curve changes from the usual 
BCS curve to a curve with strong positive curvature. 



quadratic dispersion along one direction and much flat- 
ter (higher power dependence) behavior along the other 
(transversal) direction, indicating therefore a divergence 
exponent a smaller than 1/2. In the case of saddle 
point obtained in [fj"7| , a good fit is obtained with a 
quartic dependence, leading to a sa 0.25 which agrees 
reasonably with the value extracted from the experi- 
mental H C 2 curve. We emphasize that for a given ex- 
ponent a, the normalized H C 2 curve is as universal as 
the usual BCS curve (a = 0). A suggestion of 
some sort of universality is indeed observed in Fig. 4 
of Ref. [6], where H C 2 curves of two different materi- 
als, YBa 2 (Cuo.97Zn .o3)307_5 and Tl 2 Ba 2 Cu0 6+( 5, ap- 
parently fall onto the same curve in a plot of reduced 
H C 2{T) versus reduced temperature. Such universal H C 2 
behavior is not observed in the case of a simple saddle 
point which leads to a weak logarithmic divergence in the 
density of states, pj In this case, H c2 depends on the 
coupling constant g and shows upward curvature which 
becomes stronger as g is decreased. Note that H co and 
T co have a weaker enhancement in this case, [|l4j with log 
squared relations to the inverse of the coupling constant 
g, while for the extended singularity power-law relations 
have been obtained. 

It has been assumed throughout the paper that the 
VHS was pinned at the Fermi level. It has been shown in 
Ref. juj that a deviation from the Fermi level of the VHS 
leads to H C 2 saturation at low temperature, the tempera- 
ture range of this region being proportional to the energy 
difference, kBT cross ~ Ep — E v ^. However, this has not 
been observed in the experimental i? c2 curves, |^|,[| even 
though the respective samples are in the overdoped or 
underdoped regime. According to the van-Hove scenario, 
one should expect a certain deviation of the VHS from the 
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Fermi level in these regimes in order to account for the 
decrease of the zero field critical temperature. However, 
it is possible that, (at least, in some interval of doping 
range) the reduction of critical temperature results not 
from the deviation of the VHS from the Fermi level, but, 
instead, from the weakening of the VHS due to a reduc- 
tion of the extent of the saddle point as suggested by King 
et al. [fl8| Photoemission experiments on YBa2Cu3C>6.9, 
H YBa 2 Cu 3 6 .5, and YBa 2 Cu 3 6 .3 @ support this 
picture, since they report a clear doping independence of 
the pinning of the Fermi level at the VHS. Moreover, this 
doping independence of the pinning is predicted by many 
numerical studies, from slave-boson calculations [ p4|j25| 
to renormalization group calculations prjfl . 

In conclusion, we have studied the effect of a power- 
law divergence of the density of states at the Fermi level 
N(e) ~ (T a on the upper critical field of a clean isotropic 
weak-coupling superconductor. We have shown that for 
a weak divergence (a less than 1/2), the zero tempera- 
ture critical field is finite, but strong positive curvature 
appears in H C 2 as a approaches 1/2. For a stronger di- 
vergence {a larger than 1/2), iJ c2 (T) has a power-law 
divergence at T = 0. A very good one-parameter fit was 
obtained to the experimental results by Mackenzie et al 
j|. According to the picture described in this paper, the 
anomalous H C 2 behavior reflects the short-range enhance- 
ment of the pair propagator and the unusual temperature 
dependence of the thermal length which result from the 
existence of a strong divergence of the density of states 
at the Fermi level. 
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